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ABSTRACT 

Martian craters In the size range 10 to 250 km follow 
a }og-nortna1 size-frequency distribution law more closely 
than the heretofore accepted linear relationship between the 
log of the crater diameter and the log of the crater density. 
Analysis techniques based on the log-normal model yield pos- 
sible evidence for the size-frequency evolution of crater - 
producing bodies. 

Some regions on Mars display excessive depletion of either 
large or small craters; the most likely causes of the depletion 
are considered. Apparently, eollan sedimentation has markedly 
altered the population of the small craters south of -30° lat- 
itude. The general effects of crater obliteration in the southern 
hemisphere appear to be confined to diameters of less than 20 km. 

A strong depletion of large craters in a large region just south 
of Deuteronllus Mensae, and In a small region centered at 35° 
latitude and 10° west longitude, may Indicate locations of subsur- 
face Ice. 
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BACK6R0UND OF CRATER SIZE-FREQUENCY STUDIES 

Thomas L. MacDonald (1931) was apparently the first to suggest that 
for lunar craters the log of the crater density is a linear function of the 
log of the crater diameter. Since that time virtually all crater size- 
frequency population analyses of lunar and martian craters have asswned 
the explicit validity of this log-log (LL) relationship. Most researchers 
have attributed any deviations from this presumed initial distribution of 
diameters to physical processes acting subsequent to crater formation. 
Plotrowski (1953) and J. Jones (1968) developed a posteriori theories to 
explain the Inevitability of the LL relationship for asteroid diameters, 
although the physical processes that they invoked were admittedly not 
realistic. The underlying assumption was that in an asteroid-asteroid 
collision all of the kinetic energy is utilized in the complete pulveri- 
zation of the entire smaller body plus some fraction of the larger body. 

The products of pulverization undergo no subsequent collisions, and only 
one body, the undisturbed fragment or the larger asteroid, remains to 
Interact with yet another asteroid. Considerable mathematical simplifications 
yield an asymptotically LL distribution of asteroid masses which is then 
superimposed on a planet, leading to a LL crater size-frequency distribution. 
The authors invoked this model because they could find no method to describe 
the size-frequency distribution of the fragments resulting from such col- 
lisions. 

Kolmogoroff (1941) considered the problem of the observed log-normal 
(LN) size-frequency distribution of comminution products. He made two 
assumptions about the breakage process and an asymptotically log-normal 
size-frequency distribution function resulted from his analysis. The 


assumptions are (a) In any interval of time all bodies have an equal pro- 
bability of being broken, and (b) the fractional sizes of the products are 
random and independent of the past history of the bodies. These postulates 
may not be the only ones that lead to a LN size-frequency distribution, but 
the assumptions are physically plausible for asteroid -asteroid collisions. 
Observations of the statistical distribution of large craters on Mars verify 
the LN nature of the size-frequency distribution and are therefore consistent 
with Kolmogoroff *s model (Woronow, 1975). 

THE DATA SET 

Unless specifically stated otherwise, all analyses will employ a data 
set consisting of nearly 10,000 martian craters in a latitude band from +65® 
(North) to -65® (South), and with diameters in the range 10 to 250 km. These 
craters were hand-measured on USGS quadrangle maps, refined to either the con- 
trolled or semi -control led level. All data were stored by 5® latitude x 5® 
longitude blocks on a PDP-ll/DEC tape file, and processed with a PDP-11 
conputer. For comparative purposes. Figure 1 shows the data plotted on log- 
log axes in a manner analogous to that suggested by K. L. Jones (1974); 
and it agrees well with the data set which he presents. This plot would 
commonly be interpreted as consisting of two straight line segments of the 
form 

LOG (crater density) = m • LOG (crater density) + b 
which intersect at approximately 20 km diameter: the large craters lying 
along a -3 slope (m * -3), which is presumed to represent the form of the 
initial distribution, and the small craters lying approximately along a 
-2 slope, considered an effect of erosion (e.g. Hartmann, 1973). Saturation 
in the 10-250 km diameter range is not q problem (Hartmann, 1973, Figure 3). 
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NATURE OF THE LOG-NORMAL DISTRIBUTION FUNCTION 

Undoubtedly, the wide acceptance of the LL based analyses Is due in 
part to the graphical ease and mathematical simplicity it offers. An under- 
standing of the LN distribution function's properties, origins, and uses is 
essential if the analytical techniques are to be grasped intuitively. Only 
the most salient points are sunwarized in this section; a thorough treatment 
can be found in Aitchison and Brown (1973); also Epstein (1947) presents a 
helpful review of Kolmogoroff 's paper on breakage products. 

Curve A in Figure 2 is a LN frequency curve. For any value of the 
abscissa the curve represents the frequency of occurrence of that value. The 
distribution has been normalized so that the area under the curve equals 
unity (i.e., any sampled value of the variate has a probability of one that 
its value will lie betv/een zero and plus infinity). 

For crater statistics, the LN distribution has two important differences 
from the LN distribution; (a) because the area under the LN curve can be 
normalized, the total number of craters is necessarily finite. (The LL 
relationship, in its broadest interpretation, predicts an infinity of craters), 
(b) At zero diameter the number of craters goes to zero in the LN model but 
not in the LL model. In both respects the LN distribution corresponds more 
closely to physical reality. 

Curve B in Figure 2 is a normal or gaussian curve. To transform a LN 
frequency curve to a gaussian frequency curve, one need only take the log- 
arithm of the abscissa values and plot the resultant values against their 
normalized frequencies. Following this transformation, the full pwer of 
the techniques designed for the analysis of gaussian data, including the 
F-test and t-test, are available for comparing one population to another 
or for comparing observed distributions to theoretical models. As with a 
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gausslan distribution function, the LN distribution function is completely 
specified by just two parameters, the mean and standard deviation, or alter- 
natively, the mean and variance. 

The reproductive properties of the LN distribution foHow directly 
from those of the gauss i an distribution. If two independent gauss fan variates 
are added the sum is a gaussian variate; consequently, if two independent LN 
variates are multiplied the product is an LN variate. Also, a LN variate 
can be multiplied by a constant or raised to a constant power and reproduce 
a LN variate. However, for all these cases the resulting mean and standard 
deviation may differ from that cf the initial variates* distribution function. 
The reproductive properties allow the assertion that a LN distribution of 
asteroid diameters will lead to a LN distribution of crater diameters. The 
reasoning is as follows: 

1. If the radii of the asteroids are LN, so are their volumes. 

2. If the asteroids within the size range of interest have' 
nearly equal densities, their masses will be nearly log- 
normally distributed. 

3. If the velocities of the Impacts are nearly uniform, or if 
the velocities are log-normally distributed and independently 
distributed with respect to the mass, then the product of 
powers of mass and powers of velocity will be LN. (A LN 
distribution of impact velocities is implied by the follow- 
ing traits: (a) The frequency of impacts must go to zero as 
the velocity of impact goes to zero, (b) Negative velocities 
are not allowed, (c) The velocity distribution Is likely to 
be positively skewed, I.e. have a long tail where fewer and 
fewer Impacts occur as velocity increases, once the value of 
the mode is passed). 
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4. Finally, If the crater diameters scale according to a law 
of the form 

Crater Diameter » K, (asteroid mass) 2 x 

(asteroid velocity) ^ 

(Gault, 1970) where the K's are constants over the range 
of crater diameters of Interest, one may conclude that 
the crater diameters should have a LN distribution. 

LOG-CUMULATIVE FREQUENCY PLOT 

The log-cumulative frequency plot Is a graphical convenience whereon 
LN populations appear as straight lines. (In this study the term "population" 
refers to the entire range of Information, zero to plus Infinity, while the 
term "data" refers to the Information in the range actually sampled. The term 
"transformed", preceding either "data or "population", means that the logarithm, 
base 10 in this study, has been taken). Figure 3, curve A, illustrates a 
hypothetical LN population. The abscissa is the log of the crater diameter 

IS t-TCO 

and the ordinate/such that a normal distribution will plot as a straight line. 
Standard works such as Abramowitz and Stegun (1972) give appropriate tables, 
series expansions and analytical approximations for the construction of thes^ 
graphs and the subsequent analysis of the data. 

Because the mean, median and mode all correspond in a gaussian distri- 
tHJtion, the 50 percentile intercept (median) is also the mean (y) and the mode 
of the transformed population. The slope of the line indicates the standard 
deviation of the transformed population (a), which can be calculated as follows; 
Read the values at 16%, 50% and 84%, calling them vl6, v50 and v84 respectively. 
Substitute into Equation 1. 

0 - l/2(v84 - vl6) (1) 

For the hypothetical population Illustrated in Figure 3, read the following 
values: * 


16 percentile - 1.32 
50 percentile - 1.65 
84 percentile - 1 .98 

Therefore I the transformed mean, mode and median are 1.65, and the standard 
deviation of the transformed population Is at 0.33. Equation 2a gives the 
most numerous diameter, or the mode, of the untransformed population as 25.1. 
Equation 2b gives the median of the population as 44.7. 

Mode = e (2a) 

Median = 10^ (2b) 

where A “Ml O') 

Unfof tunately, observations will always be restricted to some relatively 
narrow range of crater diameters, bounded by our Interest, Image resolution, 
planet size or population homogeneity. In such cases, the 50 percentile value 
Is only the median of the transformed data; It Is not the mean of the trans- 
formed data, nor the mean of the transformed population. The slope Is still 
an Indicator of the standard deviation of the transformed data, by Equation 1, 
though It does not indicate the standard deviation of the transformed popula- 
tion. Figure 3, curve B, shows the crater data of Figure 1 as It appears when 
plotted on log-cumulative frequency axes. 

An especially useful property of the log-cumulative frequency plot Is 
that It does not display crater density. (Crater densities and crater density 
differences are still Important parameters, but they can be analyzed Indepen- 
dently). Figure 4 has the same 10,000 craters that are shown In Figure 3, 
along with two other lines consisting of 50% and 5% of that data selected 
randcxniy. As long as the data are drav/n frexn the same population, they will 
graph nearly Identically. Only the probable errors associated wUh the mean 
and standard deviation of the transformed data change with the size of the 
data set; that Is to say that the estimators of the mean and standard devla- 
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tlon are consistent . Figure 5 Is the corresponding LL version of Figure 4. 

The same data are in both of these figures, but the crater density differences 
attract the eye In Figure 5, and the fact that the samples are all drawn from 
the same population Is difficult to see. In this artificial case, one knows 
the adjustment factors (x2 and x20) that would shift the 50^ and 5 % curves 
back to coincidence with the 100% curve. In practice though, 1f two surfaces 
had different crater densities, one would not know which curve to adjust nor 
by how much. A further benefit of analysis based on the log-cumulative fre- 
quency plots Is Immediately apparent: reconnaissance surveys of regional 
population parameters can be quickly accomplished by randomly selecting a 
relatively small portion of tfje entire crater population available and measur- 
ing and testing only those selected craters. However, one must be able to 
accept the larger uncertainties associated with the smaller sample sizes In 
order to gain a more rapid analysis. 

TESTING THE ALTERNATIVE MODELS 

Figure 2 Illustrates that a large portion of the transformed data plots 
on a straight line. This strongly suggests that a LN model 1s appropriate. 
Hwever, at small diameters the line turns dowm^rd, and off the top of this 
plot one v/ould find that the line turns upward. Both of these curved segments 
are largely due to truncation effects associated v/lth the 10 and 250 k«i diameter 
limits, although at least some of the deviation at the small dian^ters may well 
be due to erosion of craters with diameters less than 20 km. 

In order to explain the more rigorous statistical tests that verify the 
log-normality of the data, we must use a mathematical model appropriate to a 
doubly truncated log-normal distribution. 

tet X - L °3 - P 

where D Is the crater diameter, y is the transformed population mean and o Is 
the transformed population standard deviation. 
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The probability that an observation, X, will be less than or equal to 

X Is 2 

(X fx) = 1 / 2 ir| e"^ dt 

*flO 

If there were no truncation of the population. 

Let be the minimum diameter and be the maximum diameter considered, 


then 


X, . ‘■"a °1 • y 


. Log D 2 - M 
*2 


The probability distribution function for the truncated LN distribution Is then 

Pr ’(XSx) - Pr(X<Xn) 

* FCXSx^r- ' FrIXix, ) 

X rt Xi 


. f e-t'/'dt - f 

■ JmO» 

fXz *2/2 .t^/2 

e ^ dt - e ' ^dt 
J «pOft J mm 


(3) 


In order to determine If the data can be satisfied by a doubly truncated 
log-normal model, as given In Equation 3, we must first determine values of u 
and 0 which are appropriate to the transformed population. Two techniques 
have been used to accomplish this. The simplest Is devised for fitting data 
to doubly truncated normal distributions from knowledge of the first and 
second moments about the lower truncation point (Cohen, 19S0). The second 
method Is slightly more flexible and not as sensitive to non-randcxn errors; 

It Is a nonlinear regression by least squares on Equation 3. 

In either case, once the values of u and a are estimated, the test of 
Kolntjgoroff and Smirnoff (Massey, 1951; Miller and Kahn, 1962) can test the 
data against the model. Like the chi-squared test, this Is a general, non- 
parametrlc test (test results do not depend on the form of the underlying 
distribution function). The Kolirogoroft -Smirnoff method Is favored over the 


- 11 - 


chi-squared method for two reasons. First, It can provide graphical displays 
of confidence Intervals on log-cumulative frequency plots, and second, the 
method Is more powerful than the chi -squared test. 

Figures 7 through 9 compare both the LL model, with a constant slope of 
-3, and the LN model, from a least squares regression, to the actual crater data. 
The diameter ranges 10, 20, 30, and 40 km to 250 km are treated In successive 
figures. Each of the models Is represented by Its 95jS confidence bounds (I.e. 
data drawn fro ; n given size-frequency distribution model, LL or LN, will have 
only a 5X chance r violating the appropriate bounds at any point). In every 
case the data arc consistent with the LN model, but Inconsistent with the constant 
slope LL model. One might hope to salvage the LL model by c<md)1n1ng tvw> linear 
serpents, but, as shown In Figure 9, even the large diameters are not consistent 
with a line having a single slope. (A check of Figure 1 will further substan- 
tiate this). Therefore, the LL model Is unsatisfactory. (Also, two linear 
segments allow four free parameters with which to fit the data; the LN model 
requires only two; thus, by Ockham's Razor alone the LN model Is preferred). 

The fitted values of the mean and sta- Jard deviation for the LN models 
differ for each of the diameter ranges shown In Figures 7 through 9. Figure 
10 shows how the fitted parameters vary with minimum diameter. The predictions 
level out at values of u - 1.363 and a ■ 0.290. A con^arlson of the Kolmogoroff- 
Smlrnoff 95X confidence Interval appropriate to this model against the same 
diameter ranges as above, shows that good agremnent occurs for the 20 km and 
larger values of the minimum diameter, but the 10 km minimum data Is Incon- 
sistent with this model. 

POPULATION TIME EVOLUTION; GENERAL CONSIDERATIONS 

From the end of the planetary accretion until the present, the number 
of cratering events per unit time has decreased, probably more or less mono- 
tonlcally. But Is It possible that the population parameters of the Impacting 
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bodies have also changed? It can be assumed that as asteroid-asteroid collisions 
continuously occur» the mean radius of the population would continuously decrease 
(Independent of the form of the underlying size-frequency distribution). The 
standard deviation would generally Increase with time. Unfortunately, Hars 1$ 
not the Ideal locality to search for these time variations because of the action 
of size-dependent erosive and obliterative p.ocesses. Nonetheless, on Kolwo- 
5 «roff's hypothesis the development of a log-normal asteroid size-frequency 
distribution Is coupled to such a time evolution. Neukum et al. (1975) attempted 
to observe evolution effects but was unsuccessful. Because their negative result 
poses a significant threat to a LN model, I feel a critical review Is necessary. 

The data base of Neukum et al. consisted of small craters from selected lunar 
regions. However, their data covered different size ranges In different study 
regions, and to establish a "calibration curve", virtually all of their data 
required an adjustment In terms of crater density. This required averaging den- 
sities where the diameter ranges overlapped. In order to determine an adjustment 
factor. The process of averaging diminishes the possibility of recognizing 
population variations, especially on log-log axes. Their averaging schane was 
not successful; It left the crater density at any diameter In doubt by approxi- 
mately iSOrt, due to a persistent trend In their data, regardlers of size range, 
to display a shallower slope at the larger diameter extrane than at the smaller 
diameter extreme. The averaging, therefore, ccmmonly occurred between tw seg- 
ments having different slopes - leading to a large probable error in the adjust- 
ment factor. Although Neukum et al. (1975) were careful to remove suspected 
secondary craters fran their counts, the tendency of the small craters to lie 

on a slope approaching -4 may Indicate that many secondaries went unrecognized. 

In fact, their calibration curve looks strikingly like the lunar curve for 
secondaries plus primaries presented by Shoanaker (1965, Figure 47). Neither 
their data nor their technique was truly suitable to the search for tlnw de- 
pendence and they found none. 
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For the reasons mentioned previously, a Mars-based data set offers little 
> more hope of uncovering time dependence, and. Indeed, I found no Irrefutable 
cases. But finding likely regions to examine for time dependant effects or for 
abnormal eroslonal effects (treated In the next section) employs the same tech- 
nique which is introduced here. The method uses a IN model with a p and a o appro- 
priate to the planetwide population and searches out regions that violate the 
corresponding Kolmogoroff-Smirnoff bounds. To accomplish this end, each 5® x 5® 
square was tested In succession. Ccmmonly, a single square contained an Insuf- 
f1c1ent/d>f craters to warrant testing; however, when combined with enough neigh- 
boring squares, a meaningful test was possible. Therefore, the procedure employed 
requires a mlnirunt number of craters per test (minimum statistical resolution) 
and allows fre^j adjustment of the surface resolution. Figures 11 through 14 
display some of the results- of this procedure. When a violation of the limits 
occurred, the computer placed a dot In the center of the square where the analysis 
was Initiated, and drew a box around the area Included In the statistical test. 

As explained In the figure captions, one can display areas violating either the 
upper or the lower bounds over any crater diameter range and at any confidence 
level. The meandering lines In these figures (modified after Carr, et al., 1973) 
separate the plains units (predr nanily of volcanic origin) frcxn the cratered 
units. The Rcxnan Numerals used In the remainder of the paper will refer to the 
regions designated cn Figures 11 through 14. 

Having located the abnormal regions, comparison of their data by non- 
parametrlcal techniques can proceed. The method for comparing the data mean 
values Is Fisher's randomization test as described by Snedecor and Cochran (1973). 
Comparison of the data standard deviations Is with the Miller Jackknife test as 
described by Hollander and Wolf (1973). 
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Indlcatlons of possible time variations come fron Figure 10. Below 
25 km the fitted value of the mean decreases (and the standard deviations 
increase) as successively smaller crater diameters are included, demonstra- 
ting that the population parameters appropriate to those craters greater than 
25 km will not suffice upon inclusion of diameters less than 15 km. Because 
the cratered southern hemisphere greatly daninates these statistics, and there 
obliteration of the small craters presumably proceeds more rapidly than that 
of the large ones, the small craters have a lesser mean age than do the large 
craters (Hartmann, 1971). Being more recent, the small craters should represent 
a population with a smaller mean value; this is consistent with the above obser- 
vations and consistent with a time dependence of population parameters. 

Further arguments for a time evolution appear in a later section on the 
analyses of abnormal Region II (Figures 11 and 12) which shows an overabundance 
of small craters and is a relatively youthful surface. 

REGIONAL CRATER OBLITERATION 

Two general classes of crater obliterative processes mu:t be considered: 
those that remove craters with rates directly related to crater diameter and 
those with rates inversely related to crater diameter. Processes in the latter 
class are more numerous, and include most types of sedimentary filling and magmatic 
inundation. These processes undoubtedly dominate the former class in terms of 
the total effect on the crater population. Even so, their effects seem mainly 
confined to the less than 15 km diameter range, as inferred from Figure 10 which 
shows that above approximately 15 km minimum diameter, a single set of population 
parameters suffices to describe all larger diameter data sets. Furthermore, 
because a log-normal model fits even the 10-250 km data set, crater obliteration 
at the small diameters is probably not as extensive as studies using a LL model 
had indicated (e.g. K. C. Jones, 1974, Hartmann, 1973). 
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Locally, a relatively extensive obliteration of small craters will lead 
to a violation of the lower Kolmogoroff-Smirnoff bounds. Figures 13 and 14 
show just such areas; Regions IV, V, VI, VIII, IX, XII and XIII. 

Regions associated with the cratered southern terrains and with a relative 
overabundance of small craters command particular interest. If substantial ice 
exists anywhere in the martian subsurface, it may result in the more rapid iso- 
static adjustment of craters than in regions devoid of subsurface ice. Isostatic 
adjustment is probably the main process with a rate related directly to diameter 
and capable of affecting substantial surface area. In Figures 11 and 12 two such 
regions appear: Regions VII and XI. 

The other regions not yet* mentioned are Region X, the Chi*/se basin, where 
the Viking softlander is to touch down; Region II, including most of the volcanic 
plains of the northern hemisphere; and Region III, including the volcanic plains 
of Hesperia Planum. Region I is utilized as a reference surface and does not 
appear on Figures 11 through 14. 

A discussion of each of these thirteen regions now follows. The geologic 
descriptions of each region are from Carr, et a1. (1973). Topographic infor- 
mation is from Christensen (1975), and the isostatic states are from Phillips 
and Saunders (1975); these sources should be consulted for additional descrip- 
tions and details. 

REGION I 

Region I is defined to be the area south of the boundary between the 
cratered units and plains units and north of -30® latitude. This region is 
not an abnormal region, but a region against which to compare the abnormal 
regions. It consists almost uniformly of heavily cratered units, with ti»e 
exception of Hesperia Planum. In order to avoid the eolian blanketed terrains 
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Biapped by Soderblom, et a1 . (1974), only terrain north of -30° latitude was 
Included. However, Region I contains significant portions of abnormal Regions 
II, VII, VIII, IX. and XIII; but Region I*s data mean value is 1.371 and 
standard deviation is 0.244 compared with values of 1.386 and 0.244 for the 
entirety of the cratered units. The difference in the mean values is not 
significant at the 90% confidence level . 

Region I appears to be as good a touchstone for comparative analysis as 
available in the martian southern hemisphere. 

REGION II 

Region II is the largest contiguous anomalous region identified on Mars, 
covering most of the surface north of the boundary between the cratered units 
and plains units. Although this region extends south of the boundary in some 
localities, for the purposes of this study only that portion lying north of the 
boundary constitutes Region II. 

The geologic units of Region II are all lightly cratered and, therefore, 
relatively young. The units are of mixed origin, with volcanic and eolian 
deposits. The range of altitude in this region exceeds that of any other region 
studied, reaching from -2 km north of the crater Mie to +28 km atop Olympus 
Mons. Hov/ever, most of Region II lies at moderate altitudes on relatively 
shallow regional slopes. 

Region II consists of a nearly homogeneous crater population. To check 
for anomalous regions v/ithin Region II, the data from this region alone were 
fitted to a LN model, and subregions violating the Kolmogoroff- Smirnoff bounds 
for that model v/ere sought. The results appear in Figure 15. Areas around 
the volcanic cones show a significant excess of small craters compared to the 
rest of Region II. This may be attributable to many local endogenetic calderas 
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or to significant age differences. A similar test for an excess of large craters 
revealed that north of the boundary between the cratered units and plains units, 
only Region VI has a significant excess of large. craters. 

Region II, considered as a whole, derives Its anomalous nature from an 
overabundance of small craters with respect to the planetwlde average. The 
region 1$ well defined In the 10 and 25 to 250 km diameter range anomaly maps 
(Figures 11 and 12). Erosion of craters and time evolution of Impactors causes 
this persistent overabundance of small craters. The erosion Is not manifested 
in Region II; a relative lack of erosion there compared with the southern hemi- 
sphere causes the excess of small craters observed In Region II. 

S(Uiie evidence supporting the time evolution of Impacting bodies has already 
appeared earlier In this study. The examination of Region II adds two more 
lines of evidence: (1) Erosional effects, planetwide, are inappreciable above 

15 km diameter (see arguments at the beginning of this section); yet Region II 
continues to register a crater anomaly In the range 25-250 km diameter (Figure 
12). Erosion alone, therefore, will not explain the anomalous behavior of 
Region II. However, one expects this apparent overabundance of small craters 
In all diameter ranges if the mean diameter of the Impacting bodies has de- 
creased with time. (2) The Region I population parameters for the range 10-250 
km diameter are y = 1.228 and a ~ 0.335 where those for Region II are y = 0.975 
and a = 0.381. Therefore, 80% and 50% respectively of the fitted populations 
He In the diameter range 10-250 km. Therefore, to go from a Region II popu- 
lation to a Region I population requires removal of at least 40% of the small 
craters. This much (minimum) obliteration is impossible to confine to the 
10-15 km diameter range; the simpler alternative of time dependent population 
parameters Is favored. 
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Region III 

Region III essentially coincides with Hesperia Planum, and like much of 
Region II this region consists of volcanic plains. The crater anomaly maps 
in Figures 11 and 12 illustrate the close kinship of these two regions, and 
a separate discussion of Region III will not be undertaken. The observations 
and arguments pertinent to Region II also pertain to Region III without sig- 
nificant alternation. 

REGIONS IV AND V 

These two regions, IV and V, are considered together because of their 
similarity in geologic and geographic setting, crater densities, and data mean 
values and standard deviations. Both regions lie south of -30° latitude in 
the mantled terrains of the cratered units (except where Region V breaks into 
Hellas Planitia, which, from a lack of craters there, is statistically insig- 
nificant). 

A dearth of small craters up to 20 km diameter characterizes these regions 
on Figures 13 and 14. (Intermediate diameter range anomaly maps confirm the 
persistence of these regions). Although Soderblom et al. (1974) did not find 
much evidence for mantling of craters with diameters greater than 10 km, the 
effects are present. The favored interpretation of these and intermediate 
crater anomaly maps is that severe eolian sedimentation has affected all of 
Regions IV and V to beyond a diameter of 10 km, and that in some locales the 
blanketing is of importance to beyond 20 km diameter. 

Recent local disturbances caused the gaps between these two regions. The 
gap at 25° West longitude has suffered the effects of the Hellas event, which 
raised high mountains and ejected considerable debris, undoubtedly obliterating 
large and small craters alike. Similar arguments pertain to the formation of 
the Argyre basin and the gaps which surround it. 


rf-V. 
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REGION VI 

Region VI lies inmedlately to the north of the northernmost extent of 
the cratered units. It consists of plains units and knobby terrain. The 
imagery in this area is among the poorest obtained anywhere on Mars; therefore, 
conclusions concerning Region VI are necessarily tentative and correspondingly 
brief. 

The crater anomaly maps (Figures 13 and 14) for Region IV reveal an appar- 
ent depletion of craters in the diameter range 10 to 15 km. Because this region 
is in the northern mantled terrains, the small craters could have been inundated 
with dust. This conclusion remains tentative because of the poor quality of the 
available data; perhaps image resolution alone could be blamed for the ancmtaly. 
REGIONS VII AND VIII 

Astride the border between the cratered units and plains units, southeast 
of Acldalum Planitia, lie Regions VII and VIII. These two regions are discussed 
together because they hold a large portion of their areas in common. These 
regions possess no unusual geology. The area Is rather flat and In isostatic 
equilibrium. The crater density in this area, however. Is below that of Region . 
I by 9 standard deviations 1n the 10-250 km diameter range. (This is calculated 
with the portion of Region VII lying north of the boundary between the cratered 
units and plains units excluded; even larger differences result if this northern 
portion Is Included). 

The crater statistics of this area are extremely complex. Briefly, from 
Figures 11 through 14 and intermediate diameter range maps. Region VIII shows 
a relative depletion of small craters (mostly In the range 10 to 20 km diameter). 
This is a very severe depletion, judging frwn the 9 a crater density difference 
with Region I. In addition. Region VII Indicates that the northwest portion of 
Region VIII also has a relative overabundance of craters in the mid-diameter 
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ranges (20 to >25 km) with respect to the large craters. Even so, the absolute 
crater densities are low In the midrange of diameters (at the 95X confidence 
level} when compared to Region 1. 

The proposed explanation of these observations requires either crater 
obliteration at all diameters, but with greatest efficiency for the smallest 
and largest diameter extremes, or c^plete obliteration of all craters with 
subsequent obliteration of the largest and smallest diameters of the newly 
formed population. No single erosive process can meet these requ1r@nents. 
Because Sagan et al. (1973) found few wind streaks In this area, imich eollan 
sedimentation Is unlikely; an Interpretation consistent with Soderblom's et al. 
(1974) observation of a lack of mantling of the small craters there. More 
likely, most of the small craters have suffered magmatic Inundation. A few 
examples of what appear as breached crater walls (e.g. the crater Gill and the 
small crater north of It) lend support to this Interpretation. 

The agent responsible for the obliteration of the large craters Is of 
great Interest. The difference between this area's large-crater population 
and that of Region I Implies that approximately 10% of the craters greater than 
25 km have been erased. A real possibility exists that this fraction of craters 
has drowned In eruptions emanating fran their own floors or nearby sources. 

But the proximity of this area to the knobby terrain (Immediately to the north- 
east) fuels an alternative explanation: the large craters have undergone rapid 
Isostatic rebound due to the presence of subsurface Ice. Both Sharp (1973) 
and Milton (1973) argue that the knobby terrain may have resulted from the 
presence of subsurface Ice. 

The coupling of an active volcanic province with subsurface Ice Is not 
Impossible. The transition of this area from a cold to a thermally active 
province could have caused the collapse Into knobby terrain and the very rapid 
Isostatic rebound of large craters. 
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REGION IX 

Region IX lies just east of Region VIII and entirely within cratered 
units. It rests on a north-facing slope spanning a 2 to 5 km altitude range. 

Figures 12 and 13 show Region IX to have a marked deficiency of small 
craters. This trend is more firmly established in Region IX than in any other 
region found on Mars, maintaining the deficiency to beyond 25 Ion dleuneter. 

The cause for such a pronounced anomaly 1s unknown but the proximity of this 
region to Region VIII strongly suggests that magmatic Inundation could be the 
active agent here as well. 

REGION X 

Region X lies adjacent to, though not conterminous with. Region II. It 
Is predominantly north of the boundary between the cratered units and plains 
units and contains the basin of Chryse Planitla. As In Region II, Region X 
consists primarily of plains un1ts» although It also Includes some chaotic and 
cratered terrains In its southern and eastern extents. Topographically, Region 
X lies mostly within a basin which descends to 3 km below the adjacent cratered 
terrains. The Chr/se basin exhibits a negative Isostatic anomaly suggesting the 
presence of a thick, low density unit. Channels connect the basin with the 
adjacent chaotic terrains, although flow direction indicators are ambiguous. 

Region X differs from Its neighbor. Region II, In at least one major statls 
tical feature: In Region X the anonaly appears only on the 10-250 km ananaly 
map (Figure 11), It Is not present on a 15-250 km anomaly map, whereas Region 
Il's anOTialistic behavior extends to the large diameter ranges. However, a 
complete absence of large craters in the basin has caused a considerable loss 
of ground resolution, and thus the adjacent normal terrains were automatically 
Incorporated In these statistical tests. This incorporation obscures any 
anomalistic population within the basin. . 
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The similarity of the crater populations of Regions X and II requires 
only that both regions be erased of all craters at approximately the same time 
and then be left to start accumulating craters anew. Thus, the formation of 
Oiryse basin must be contemporaneous with or predate the general obliteration 
episode which affected the entire northern plains. Two plausible processes 
could have delivered sediments to the basin, obliterated the craters, and 
caused the negative anomaly: (1) eollan transport from the North, and (2) 
aqueous transport fr«n the South and West. 

Although Soderblom et al. (1974) did not Include Chryse basin within their 
mantled terrains, Sagan et al. (1973) mapped strong wind direction Indicators 
coming from the Northeast, off Acldallum Planitia and the mantled terrains. 

Thus at least some craters may have suffered obliteration by eollan deposition. 

The existence of aqueous agents Is more problematical. If the chaotic 
terrains have resulted from the collapse of overlying sediments after the with- 
drawal of subsurface Ice, then any liquid produced would have come down Into 
the basin (assuming conditions were correct for melting rather than for subli- 
mation of the Ice). This could explain the channels connecting the basin and 
the chaotic terrains. 

Whether one or both of these plausible eroslonal agents actually occurred, 
the crater population in Region X Indicates that a major episode of obliteration 
there and in Region II were roughly coincident In time. 

REGION XI 

Region XI is a small region lying northeast of Argyre basin In cratered 
terrains of normal appearance. It rests on a shallow slope facing the basin 
and elevated above the basin floor by approximately 3 km. 

Although Region XI's crater size-frequency distribution Is interesting 
fran the standpoint of subsurface ice. Its small areal extent poses major prob- 
lems. The crater anomaly there may simply be a coincidental grouping of craters 
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even If It is not* the data set and surface area are small and not readily 
amenable to detailed analysis. 

At tKeSOX confidence level no differences in the data mean values or 
standard deviations for the size range 10-25 km exist between Region XI and 
I. In the diameter range 25-250 km. however, the mean value of Region XI Is 
less than that of Region I at the 99X confidence level. This Is consistent 
with abnormally rapid isostatic rebound and concomitant obliteration of the 
largest craters. Unfortunately, the region is spar.e in morphologies generally 
considered as indicators of possible sites of subsurface ice. Thus, the likeli- 
hood of subsurface ice in Region XI remains poorly assessed. 

REGIONS XII AND XIII 

Regions XII and XIII are presented only for the sake of completeness. Both 
of these regions seem to be statistical quirks rather than meaningful anomalies. 
These regions show a depletion of craters in the midrange diameters (20 to 30 km) 
only. A rather complex geological history would need to be constructed, with 
superpositioning of two or more obliterative processes, in order to explain the 
observed anomalies. With the requisite superpositioning of processes, at least 
one process would be expected to have also affected the craters adjoining these 
regions. Therefore, regions with different types of anomalies would border 
Regions XII and XIII. Because no such “paired** anomalies are found, the favored 
explanation is that of statistical coincidence. 

CONCLUSIONS 

To analyze a crater population that results from an integration of a time 
dependent Impacting population is a difficult task. Added C(xnpli cations arise 
because different diameter ranges may be under the dominion of different physical 
processes and are, therefore, likely to possess different distribution parameters 
or even different distribution functions. By restricting the range of crater 
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diameters and by considering only the extreme wing of the size-frequency 
distribution, some gains have been made In obtaining the form of the distribu- 
tion law. A linear relationship between the log of the crater diameters and 
the log of the crater densities is a less satisfactory representation than Is 
a log-normal law, as judged by formal statistical Inference. 

Twelve regions differ significantly frotn the planetwide crater population 
average. Simple processes of sedimentation suffice to explain most of the ob- 
servations, but some evidence for subsurface Ice, and time evolution of the 
Impacting bodies also exists. 

Replacing a LL law by a LN law calls into question a great portion of 
the previously held conclusions. The Impact Is greatest on those studies report- 
ing absolute dates derived from the density of small craters. The parallel, 
straight isochrons of the LL model, must be replaced by arching nonparallel lines. 
Dating based on smaller craters may have systematically underestimated ages by 
more than an order of magnitude. The absolute ages were already poorly determined 
because of the difficulty of estimating impact rates* but the errors caused by 
application of an incorrect model are In addition to any previously acknowledged 
uncertainties. 

The relation between a bre^ .age process and the observed crater size-frequency 
distribution bears on the origin and nature of the impacting bodies. Two possible 
populations present themselves: early stray bodies In nearly circular orbits 
and asteroids In more eccentric orbits. Both of these populations probably 
once contained a sufficient density of bodies to facilitate abundant mutual col- 
lisions. Unfortunately, these two possible populations each carry separate 
forecasts for the success of Interplanetary time correlation. If cratering 
was primarily caused by bodies in low eccentricity orbits, then neither the 
crater density nor the population parameters must correlate between planets*. 
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However, If the cratering was primarily caused by bodies In high eccentricity 
orbits, then the population parameters, and perhaps the crater densities, 
should correlate between planets. 

The analyses of this study do not exhaust the possibilities of crater 
population studies on Mars, but they do Illustrate the use of new and more 
rigorous analytical techniques applicable not only to Mars, but to all cratered 
bodies. 
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FIGURE CAPTIONS 

The crater data (10-250 km diameter) used in this paper as tK«y 
appears when plotted on log-log axes in a manner analogous to 
that suggested by K. L. Jones (1974). In comparison to the line 
of -3 slope, the data are gently, but continuously, bending. 

The log-normal frequency curve (A) can be transformed to a gaussian 
frequency curve (B) by plotting the log of the diameter against 
the normalized frequency. The positive skewness of the log-normal 
curve results in the successive rightward shift of the median and 
mean away from the mode. 

The crater data (10-250 km diameter) as in Figure 1, but plotted 
on log-cumulative frequency axes. The long linear segment suggests 
the applicability of a log-normal size-frequency distribution model 
for large martian craters. A numerical example in the text utilizes 
the hypothetical population lines. 

The three curves represent 100%, or all 10,000 craters, and 50% and 
5% selected at random from the total data set. Because crater 
densities are not displayed on a log -cumulative frequency plot, the 
fact that all three data sets were drawn from the sane population 
is immediately apparent. These data sets are repeated on log-log 
axes in Figure 5. 

The three curves represent 100%, or all 10,000 craters, and 50% and 5% 
selected at random from the total. Crater density differences are 
apparent, but the fact that all three data sets are drawn fYorn the 
same population is difficult to recognize. These same three data 
sets are on log-cumulative frequency axes in Figure 4. 
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The Kolmogoroff-Smirnoff 95% confidence limits appropriate to 
a log-log distribution model (v/ith a constant ^ope of -3) and a 
log-normal distribution model (with y ■ 1.228 and o “ 0.335 as 
determined by a least squares fit to the 10-250 km diameter data) 
compared to the crater data over the range 10-250 km diameter. 

The LN model provides a good fit to the data but the LL model 
does not. 

The Kolmogoroff-Smirnoff 95% confidence limits appropriate to a 
log-log distribution model (with a constant slope of -3) and a 
log-normal distribution model (with y = 1.317 and a ~ 0.504 as 
determined by a least squares fit to the 20-250 km diameter data) 
compared to the crater data over the range 20-250 km diameter. 

The LN model provides a good fit to the data but the LL model 
does not. 

The Kolmogoroff-Smirnoff 95% confidence limits appropriate to a 
log-log distribution model (with a constant slope of -3) and a 
log-normal distribution model (with u = 1.365 and o = 0.289 as 
determined frwri a least squares fit to the 30-250 km diameter data) 
compared to the data over the range 30-250 km diameter. The LN 
model provides a good fit to the data but the LL model does not 
(note violation of the upper portion of the data curve). 

The Kolmogoroff-Smirnoff 95% confidence limits appropriate to a 
log-log distribution model (with constant slope of -3) and a 
log-normal distribution model (with P = 1.309 and o = 0.304 as 
determined by a least squares fit to the 40-250 km diameter data) 
cOTipared to the data over the range 40-250 km diameter. The LN 

t 

•.i(kI?1 provides a good fit to the data, but the LLmodel does noti 
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Figure 10: The population parameters > mean and standard deviation, 

obtained by fitting the data to a log-normal model, graphed as 
a function of minimum diameter (holding the maximum diameter 
fixed at 250 km). The two methods used in the fitting (des- 
cribed In text) agree very well above approximately 20 km minimum 
diameter at which point the curves also nearly level out. 

Figure 11: Regions on Mars found to violate the upper bounds of a 90% 

Kolmogoroff -Smi rnoff confidence Interval appropriate to a log- 
normal model with p = 1.228 and a = 0.335 when craters In the 
range 10-250 km diameter are examined (10 or more craters per 
analysis). Vastitas Borealis and Hesperia Planum have an excess 
of small craters . 

Figure 12: Regions on Mars found to violate the upper bounds of a 95% 

Kolmogoroff -Smi rnoff confidence interval appropriate to a log- 
normal model with p = 1.363 and a = 0.290 when craters in the 
range 25-250 km diameter are examined (10 or more craters per 
analysis). Several regions show an excess of small cratef^. 

Figure 13: Regions on Mars found to violate the lower bounds of a 90% 
Kolmogoroff-Smi rnoff confidence Interval appropriate to a 
’og-normal model with p = 1.228 and a = 0.335 when craters In 
the range 10-250 km diameter are examined (10 or more craters 
per analysis). Several regions show an Insufficiency of small 
craters . 

Figure 14: Regions on Mars found to violate the lov/er bounds of a 95% 

Kolmogoroff-Smi rnoff confidence interval appropriate to a log- 
normal model with p » 1.316 and a = 0.305 when craters In the 
range 20-250 km diameter are examined (10 or more craters per 
analysis). Several regions show an insufficiency of small 


craters . 
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Flgure 15: Regions on Mars found to violate the upper bounds of a 95X 
Kolmogorof f-SmI rnof f confidence interval appropriate for a 
log-normal model with y *■ 0.957 and o * 0.381 when craters 
in the diameter range 10-250 km are examined (8 or more craters 
per analysis). The values of y and o were obtained by a least- 
squares fit to the data of Regi».n II. Areas surrounding the 
large volcanic cones show an excess of small craters. 
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